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I Abstract The foreshock region of a CME shock front, where shock accelerated 

£S) i electrons form a beam population in the otherwise quiescent plasma is gcncr- 

£H | ally assumed to be the source region of type II radio bursts. Nonlinear wave 

interaction of electrostatic waves excited by the beamed electrons are the prime 
candidates for the radio waves' emission. 

■ To address the question whether a single, or two counterpropagating beam 
C^l \ populations are a requirement for this process, we have conducted 2.5D particle 

in cell simulations using the fully relativistic ACRONYM code. 
| Results show indications of three wave interaction leading to electromagnetic 

■ emission at the fundamental and harmonic frequency for the two-beam case. For 
• | the single-beam case, no such signatures were detectable. 

Qh. Keywords: Radio Bursts, Type II; Radio Bursts, Theory; Plasma Physics 

6 

c/2 ' 1. Introduction 

Solar radio bursts have been observed since the earliest days of astronomical 
radio observations. Based on spec tral morphology and d evelop ment timescales, 



radio observations, cased on spec tral m orpnoiogy and development timescales, 
| they are classified into 5 types ( Wild and McCreadv . 19501) . While suitable 



models for the emission of most types of radio bursts have since been found, 
' an explanation of type II, and to some extent type III radio bursts' emission 

I/"") ■ mechanisms remains elusive. 

\^ ' The morphology of type II bursts typically shows a two-band emission spec- 

f^) . trum, consisting of the fundamental emission band and of the harmonic emission 

' band at about twice the frequency of the fundamental (in a few, near-limb event s 

signals of third harmonic emission are also discernable, see Zlotnik et al . 1998h . 



The fundamental frequency is believed to correspond to the plasma frequency of 
the emission region, slowly decreasing over time as the coronal/interplanetary 



. tnc emission region, siowiy decreasing over time as tnc coronal/ interplanetary 

| shock travels outwards into the heliosphere. ( Cane. Sheelev. and Howard! 1987 

C$ '■ iNelson and MelroseL fl985h 
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Figure 1. Exemplary classification of solar radio bursts by their morphology in a dynamic 
spectrum. Type II bursts are characterized by their two-banded emission of fundamental and 
harmonic frequency, which is slowly drifting towards lower frequencies. 



Statistical studies have provided firm evidence of a correlation between type 
II bursts and the propagat i on of shock fronts o utwards from the sun 
( Nelson and Melrosei 1985 ; Forbes et al , 20061 ). These are in many cases driven 
by coronal mass ejection s (CME), but have also been observed as blast waves 
from flare energy release ( H. Aurass. B. Vrsnak. and G. Mann , 2002f ) . It is there- 
fore reasonable to assume that electron processes at the shock are responsible 
for the emission ( Mannl . 1995 ). At the time of radio emission, a CME shock 
front can span from the low solar chromosphere outwards to multiple solar 
radii, thus encompa ssing a large number of density scales and hence, potential 
plasma frequencies ( Pomoell. Vainio. and Kissmann . 20081) . The observed radio 
emissions, on the other hand, are quite narrowband. Consequently, the emission 
region can not be spread out along a large area of the shock front, but is rather 
required to have a small spatial scale, confined to a small subset of the shock 
(jSchmidt and Gopalswamvl 120081 ). The emission region is thus expected to have 



special properties setting it apart from the rest of the shock. 
1.1. Type III Bursts 

Very similarly to type II bursts, the so-called type III bursts also feature multi- 
banded emission, albeit with a much faster drift towards low frequencies and 
hence shorter total burst duration. 

For these bursts, conclusive observational evidence exists for a connection 
to impulsive electron release in solar flares. As the accelerated electrons travel 
outwards on open field lines, they form an electron beam population within the 
solar wind p lasma, which can exc i te Langmuir waves th rough Cherenkov-typc 
instabilities ( Karlickv and Vandas . 2007 ; Tsiklauri . 2010l ). 

These Langmuir waves (L) are t hen assumed to take part in three wave 
interaction processes ( Melrosei 19861) . leading to emission of transverse electro- 
magnetic waves (T) and sound waves (S): 



L -> L' + S 
L S + T(uj pe ) 
L' -> T(2uj pe ) 



(1) 

(2) 
(3) 
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Figure 2. Sketch of three wave interaction processes believed to be the origin of fundamental 
(left) and harmonic (right) emission of type II radio bursts. Decay of beam-driven longitudinal 
waves into sound-waves and electromagnetic waves at the plasma frequency akin to the normal 
plasma emission process is responsible for the fundamental emission, whereas a coalesce process 
of two counter-directed modes leads to harmonic emission. 



In a more detailed analysis, IWilles and Cairnsl (|2000l) treated plasma wave 
behaviour in the presence of a beam population, and identified not the Langmuir 
wave, but the so-called beam driven mode as a possible source of three-wave 
interactions. As the beams' intensity decreases, energy in this mode would be 
transferred into the background plasma's eigenmodes, thus leading to similar 
processes as outlined in eq. HEl illustrated in figure [H 

Experimental proof o f these processes has been challenging: satellite data 



LP 

(|Pulupa and Bald . I2008I ). being intrinsically point- wise, can not yield informa- 



tion about the momenta and energies of participating waves in three-waves 
interactions. 

While numerical investigati on of three- wave interaction ha s been successfully 
performed using MHD codes ijSpanier and Wisniewskil . I2011I) , that technique is 
not applicable to the kinetic waves participating in the interactions outlined 
here. Particle-in-cell codes on the other hand, which are the tool of choice for 
kinetic simulations, have been too numerically expensive to be able to treat this 
problem self-consistently in the past. With the exponential growth of computing 
power, the length- and timescales of these emission models are coming into the 
range of top-end supercomputers now. 



2. Emission Model 

Due to the similarity to type III bursts in spectral shape, it seems reasonable to 
assume that electron beams arc likewise responsible for type II radio emission. 
Assuming that the electron acceleration on the shock front is localized to a 
specific magnetic field configuration also elegantly solves the problem of the 

emissions' narrow-bandedness. 

In the predominant model, iHolman and Pessesl (|l983l) proposed that electron 



shock drift acceleration on a curved shock front is the source of these electron 
beams. Following the interplanetary field from the point where the shock is ex- 
actly perpendicular into the foreshock plasma, these electron beams then excite 
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Figure 3. Sketch of the large scale model of the proposed emission process of type II radio 
bursts. Electrons accelerated through shock drift acceleration form a beam population in the 
forcshock plasma, exciting Langmuir waves which then undergo three-wave interactions. A 
variant of this model allows for multiple sites of electron acceleration on a deformed shock 
(such as a shock ripple) injecting countcrstrcaming electron beams into the same region. 



Langmuir waves and undergo three- wave interactions much like eq. [T][3] predicts 
for type III bursts. Since the efficiency of drift acceleration increases sharply 
around the sites where the shock angle is close to 90°, the source region of the 
electron beams is confined to a small area on the curved shock, thus explaining 
the narrowband nature of the emission. 

An open question in these models remains whether the process of equation [T] 
itself - that is the decay of a forward-directed Langmuir wave into a forward- 
directed sound wave and a backward-directed Langmuir wave - is sufficient to 
produce a large seed population of countcrpropagating Langmuir waves for the 
process of equation [3] to produce significant intensities of radio emission at twice 
the plasma frequency. Alternatively, counterstreaming electron beams, directly 
exciting countcrpropagating Langmuir waves may be required, as depicted in 
figure [3] These might be created in the case of shock ripples, or deformations of 
a shock due to shock-shock interaction, as reported multiple times in cases of 
very strong type II bursts. 



3. Numerical Simulation 



In order to decide whether a two beam scenario is required for harmonic emission, 
or whether a single-beam scenario is sufficient, we have conducted simulations 
using the ACRONYM particle in cell code. 

3.1. ACRONYM Code 



The ACRONYM code, developed and maintained at the department of As- 
tronomy, University of Wurzburg, is a fully relativistic, 2nd order particle-in- 
cell code for astronomical, heliospheric and laboratory plasmas. Using MPI- 
parallelization, the code is running o n all major supercomputer platforms 
(jKilian. Burkart. and Spanierl . 120121) . 

For the simulations in this project, computing time on Jugene (BlueGene/P) 
at the Julich Supcrcomputing Centre and Louhi (Cray XT- 5) at the CSC Helsinki 
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Figure 4. Electron velocity (v x ) distribution at the start of a simulation in the single-beam 
case. A second, symmetric beam is included in the two-beam case 



was employed. The runs took a total of about 1 
machines. 



Mio. CPU-hours on these 



3.2. Setup 



Since a particle in cell code carries the intrinsic constraint of resolving the plasma 
debye length (in the case of solar wind, \d « 1cm), the maximum size of a 
simulation domain is severely limited to the resolution of plasma-microphysical 
phenomena. A simulation of the complete shock environment depicted in figurc[3] 
would be prohibitively expensive in terms of computing power. The simulations 
rather focus on the beam-plasma interaction and resulting three- wave processes 
in the foreshock plasma. 

The simulations are set up as 2.5D (2D3V) rectangular grids with peri- 
odic boundaries, which are homogeneously filled with the background foreshock 
plasma under quiescent solar wind conditions (T sa 0.5 MK, p — 2.5 ■ 10 7 cm~ 3 , 
B = 1G), with a thermal particle distribution. On top of that, either a single or 
two counterstreaming beamed electron populations are added at Wseam ~ 5^th, 
whose density is about 10% of the total electron density (see fig. HJ. The beam 
density and velocity is derived from typical Type II emission region parameters 
given in Knock et all ( 2001 ). their pitch angle distribution is centered around 
45°, following iKarlickv and Vandasl (j2007t ). 

In the two-beam case, the simulation will be naturally charge-neutral and have 
a zero total current. In the single-beam case however, the background electron 
population is suitably boosted in the opposite direction to create an appropriate 
return current, and thus also zero out the total current within the simulation 
box. 
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Table 1. Simlation parameters for the particle in cell simulation 



Physical parameters 



Background plasma density 


P 


2.5 ■ 10 7 cm- :i 


Temperature 


T 


480280 K 


Background magnetic field 


B x 


1G 


Plasma frequency 


U)p 


2 ■ 10 s Rad/s 


Beam density 


I'D 


0.1 * p 


Numerical parameters 


Box size in cells 


N x 


8192 




Ny 


4096 


Cell size 


Ax 


0.95 cm 


Timcstep 


Ar 


1.2 ■ 10~ n s 


Particles per Cell 


n 


200 Background + 20 Beam 



Energy content, two-beam 



Energy content, single-beam 
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Figure 5. Energy distribution within two-beam (left) and single-beam (right) simulations, 
over time. After about 20 plasma timescales (uip 1 ) , the beam instability causes an increase in 
electric field energy. Subsequently, the transverse magnetic field energy also rises. (Background 
B-Field intensity is not shown, as it is constant throughout the simulation) 



4. Results 



Since the simulation is started with a quiescent thermal background plasma, 
into which the electron beam populations are injected, it takes a while (about 
t = 20a;- 1 ) for the beam-driven instability to create a suitable amount of wave 
intensity. This can be easily observed by looking at the energy distribution graphs 
of single- and two-beam simulations, which are shown in figure [5j 

In this graph, total particle kinetic energy, electric field energy and the energy 
content of the transverse magnetic field components is displayed over simulation 
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Figure 6. Dispersion plot of the longitudinal electric field component along the beam direction 
in a one-beam simulation. A strong excitation of the beam-driven mode, around resonance with 
the plasma frequency is visible in both forward direction (left). 



time (the background magnetic field component is left out, since it is constant 
in good approximation). 

The total particle kinetic energy contains a large contribution from the back- 
ground plasma's thermal motion, making up about 61% of its value at the start of 
the simulation. The remaining 38%, or 4.3 • 10 4 erg (calculated from the particles' 
initial velocities) are contribution of the beam kinetic energy. 

The electric field energy reaches a maximum at t « 50o;~ , with a coinciding 
rise in energy in the transverse magnetic fields, implying an energy transfer from 
electrostatic to electromagne tic modes. 

Confirming the results of iKarlickv and Vandasl (120071) . a peak electric field 
energy of about 2.8% of the beam kinetic energy is observed in both simulations. 
The energy content in the transverse magnetic fields at the end of the simulation 
runs lies about 3 orders of magnitude below that of the electric fields. 

4.1. Single-Beam Scenario 

In order to obtain quantitative information about the intensities of individual 
wave modes within the simulation, electric and magnetic field quantities were 
Fourier-transformed in spatial and temporal dimensions. The results arc plotted 
as k vs. oj intensity plots, in which each wave mode can be identified via its 
characteristic dispersion relation. 

Figure [5] shows this dispersion plot for the longitudinal electric field com- 
ponent (E x ) along the beam direction, with k \\ B . As expected, electrostatic 
longitudinal wave modes are discernible in this plot. In addition to the Langmuir 
wave, which is forming a parabolic shape around the plasma frequency at k = 0, 
an additional, strong area of intensity is visible at negative k values. The dis- 
persion relation of this wave corresponds to a beam driven mode of the electron 
beam, which excites oscilla tions in the electron bea m's density near resonance 
with the plasma frequency ( Willes and Cairns . 2000h . 
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Figure 7. Dispersion plot of a transverse magnetic field component in the single-beam simula- 
tion, with k perpendicular to the beam direction. The electromagnetic mode is clearly visible 
as a parabola, with a cutoff at the plasma frequency (splitting into L- and R-Mode is not 
sufficiently numerically resolved). No additional unexpected wave modes are observed. 



By taking the sum of the energy content in the forward-propagating half space 
and comparing with the with the backward-propagating one, the beam driven 
modes' energy can be determined to Ef = 475 erg, whereas the backward direc- 
tion only contains Ej, = 0.3 erg of energy. Note however that the energy content 
determined from this dispersion plot only contains modes which are exactly axis- 
parallel (up to numerical precision) and is therefore grossly underestimating the 
beam-modes complete energy contribution (compare figure [5]). 

In figure [3 a similar dispersion plot for a transverse magnetic field component 
is shown. Apart from some thermal background excitation of the electromagnetic 
mode, no unexpected features are apparent in this plot. Especially, no indications 
of three-wave interaction are visible here. 

4.2. Two-Beam Scenario 

The longitudinal electric field dispersion plot for the two-beam case (figure 
[SJ shows a symmetric structure of two countcrpropagating beam-driven mode 
peaks, similar to figure [5] Additionally, a somewhat weaker feature is visi- 
ble around 2lo p , which may be indicative of the three- wave coupling processes 
outlined in section ri.il 

The dispersion plots of transverse magnetic field components (figure [5]) show 
the expected dispersion relation of the electromagnetic mode, as predicted by 
linear theory. Additional bands of intensity at the fundamental and harmonic 
emission frequency are also visible. Consistent with the predicted processes out- 
lined in section II. 11 these bands are likely created by nonlinear interaction of 
beam-driven waves with opposite direction of propagation. 

For most fc-valucs, the nonlinear product wave does not match an eigen- 
mode of the background plasma and is hence unable to excite a propagating 
wave, explaining the low intensity of the feature in the dispersion plot. For de- 
values matching the electromagnetic mode however, electromagnetic waves can 
be resonantly excited, and travel through the plasma. 
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Figure 8. Dispersion plot of the longitudinal electric field component along the beam direction 
in a two-beam simulation. A strong excitation of the beam-driven mode, around resonance with 
the plasma frequency is visible in both forward and backward direction. Additionally, weaker 
but significant excitation appears to be present around 2u) p 
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Figure 9. Dispersion plot of a transverse magnetic field component in the two-beam sim- 
ulation, with k perpendicular to the beam direction. The electromagnetic mode is clearly 
visible as a parabola, with a cutoff at the plasma frequency (splitting into L- and R-Modc 
is not sufficiently numerically resolved). In addition, two broad features at the fundamental 
and harmonic frequency are visible, which are in resonance with the em-Mode. (Intensities in 
arbitrary units) 

5. Conclusion 

We have compared particle in cell simulations of two different microphysical 
emission scenarios of type II radio bursts in a CME foreshock settings: one 
with a single electron beam accelerated by the shock, and another with two 
counterstreaming shock populations of equal strength. 

While the single-beam scenario has shown production of electrostatic waves 
through beam-driven instabilities which are in line with theoretical predictions, 
emission of electromagnetic waves with fundamental and harmonic frequencies 
which would be the fingerprint of a type II burst, have not been observed. 

In the two-beam simulation, electrostatic modes are likewise excited, leading 
to countcrpropagating wave populations. Here, indications of nonlinear wave in- 
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teraction are visible in the dispersion diagrams of transverse magnetic field com- 
ponent, with maxima at the plasma frequency and its first harmonic, conforming 
to the theoretically predicted emission process of type II bursts. 

The results presented here should not be taken as a literary proof that a 
multi-beam environment is a necessity for type II emission: The canonical pro- 
cess to create counterpropagating wave populations for subsequent interaction 
is scattering of Langmuir waves on sound waves (eq. [I]). By using a collisionless 
particle in cell approach for our simulations, we have completely neglected the 
effects of sound waves on the plasma, thus likely creating somewhat unrealistic 
wave dynamics. 

Furthermore, due to the large numerical demands of the particle in cell sim- 
ulation, the total physical simulation time was limited to a few microseconds. It 
is therefore possible that the electrostatic wave population has not relaxed into 
an equilibrium condition by the end of the run, such that scattering effects may 
have been underrepresented. 

Still, the demonstrated strong dependence of harmonic radio emissions on 
the existence of counterpropagating electrostatic waves confirms long-standing 
theoretical predictions. In plasma environments where this mechanism is of 
importance, typical damping time of these waves has to be much larger than 
scattering times if a single beam mechanism is to be responsible for the radio 
emissions. 

5.1. Outlook and further work 

In this work, only the two extreme cases of a single beam and counterpropagating 
of equal strength beams have been investigated. In an actual type II emission 
region, beams will likely have different strengths and pitch angle distributions 
if they are originating from separate acceleration regions. The next step will 
therefore be simulations with varying beam intensities and angles. 

The parameter space of these simulations further consists of the density / 
plasma frequency of the ambient plasma, magnetic field strength and temper- 
ature. Given sufficient computing time, parameters studies in this space may 
follow suit. 
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